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I.  INTRODUCTION 

Disturbances  are  defined  as  the  uncontrollable  inputs  which  act  on  a 
dynamical  system.  There  are  many  varieties  of  disturbance  inputs  which  can  be 
associated  with  a  controlled  system  and  they  are,  for  the  most  part,  com¬ 
pletely  unpredictable  in  magnitude  and  in  their  arrival  times. 

In  practice,  additive  disturbances,  i.e.,  disturbances  which  are  repre¬ 
sented  by  terms  added  to  the  plant  state  equation,  can  arise  from  motivating 
effects  external  to  the  plant  (external  disturbances)  or  from  motivating 
effects  arising  from  the  physical  characteristics  of  plant  subsystems  or 
Internal  plant  dynamics  (Internal  disturbances).  Further,  these  disturbances 
can  be  divided  into  two  categories:  (a)  noise  disturbances  -  characterized  by 
random  and  erratic  behavior  with  relatively  high-frequency  content  and  (b) 
waveform  structured  disturbances  -  characterized  by  a  degree  of  waveform  regu¬ 
larity  which  can  be  described,  piecewise  in  time,  by  differential  equations 
forced  by  sparse  sequences  of  Impulses.  The  nature  of  these  disturbances  may 
be  either  completedy  known  (through  direct  prior  or  realtime  observation  or 
test) ,  completely  unknown  (random-like),  or  partially  known. 

Johnson  (1-6,  10]  Introduced  the  idea  of  mathematically  describing 
uncertain  waveform-structured  disturbances  by  representing  them  as  a  weighted 
linear  combination  of  known  basis  functions  of  the  form 

n 

w(t)  -  Z  c±f 4( t) ,  (1) 

i-1 

where  w(t)  is  the  plant  disturbance  vector  and  is  a  p-vector  and  the 
weighting  coefficients  c^  are  completely  unknown  constants  which  can  change  in 
magnitude  in  a  random,  once-in-a-while ,  fashion.  The  basis  functions  f^(t) 
are  completely  known  because  they  are  chosen  by  the  designer  based  on  the 
waveform  patterns  exhibited  (or  thought  to  be  exhibited)  by  the  disturbance. 

Johnson  [1-11]  developed  a  control  engineering  design  technique, 
referred  to  as  Disturbance  Accommodation,  wherein  a  combination  of  waveform¬ 
mode  disturbance  modeling  and  state-variable  control  methods  are  utilized  to 
design  controllers  which  will:  (1)  absorb  (counteract),  (2)  minimize  or  (3) 
constructively  utilize  the  effects  of  uncertain  disturbances  on  the  plant. 
Three  main  classes  of  controllers  are  considered  within  the  overall  cognomen 
of  Disturbance  Accommodating  Control  Theory.  These  are,  (1)  Disturbance 
Absorption  Controllers  (DAC),  (2)  Disturbance  Minimization  Controllers  (DMC), 
and  (3)  Disturbance  Utilizing  Controllers  (DUC).  Each  class  of  controller  has 
its  own  associated  design  goals  and  design  methodology.  The  mathematical 
theories  of  DAC  and  DUC  were  thoroughly  developed  in  [1-12].  The  theory  and 
techniques  associated  with  DMC  were  compiled  and  extended  in  [13]. 

A  number  of  examples  were  presented  in  [13]  to  illustrate  the  applica¬ 
tion  of  various  disturbance  minimization  techniques.  This  report  will  present 
further  results  on  the  application  of  disturbance  minimizing  control  design 
techniques  to  linear,  time-invariant  state  set-point  regulators. 


II 


LINEAR  DYNAMICAL  SYSTEMS 


The  class  of  systems  to  be  considered  in  this  report  are  " linear,  time- 
invariant,  dynamical  systems1*,  so-called  because  the  vector  differential 
equation  for  the  state  x(t)  is  a  linear  differential  equation,  the  transfor¬ 
mation  between  the  state  space  and  output  space  is  linear,  and  the  elements  of 
the  matrices  in  the  plant  model  are  constant  with  respect  to  time. 

These  systems  will  be  represented  by  equations  of  the  general  form 

x(t)  -  Ax(t)  +  Bu( t)  +  Fw(t)  (2) 

y(t)  -  Cx(t)  +  Eu( t)  +  Gw(t) ,  (3) 

where  x(t)  is  the  plant  state  vector  and  is  an  n-vector,  u(t)  is  the  plant 
control  input  vector  and  is  an  r-vector,  w(t)  is  the  plant  disturbance  vector 
and  is  a  p-vector,  y(t)  is  the  plant  output  vector  and  is  an  m-vector  and  A, 

B,  F,  C,  E,  G  are  appropriate  size,  known  matrices  with  time-invariant  ele¬ 
ments.  In  addition,  the  general  form  of  the  disturbance  state  model  is  [10] 

w(t)  -  Hz(t)  +  Lx( t)  (4) 

z(t)  -  Dz( t)  +  Mx(t)  +  o(t),  (5) 

where  z(t)  is  the  p-  dimensional  disturbance  state  vector,  o(t)  is  a  sparsely 
populated  vector  impulse  sequence  and  H,  L,  D,  M  are  appropriate  size,  known 
matrices. 


III.  BACKGROUND 


In  [13]  several  methods  were  presented  for  minimizing,  via  direct 
control  action,  the  effects  of  constant  disturbance  components,  which  are  not 
completely  absorbable,  on  linear,  time-invariant  state  set-point  regulators. 
The  metric  used  for  the  minimization  process  is  the  norm  defined  by 


2 

Q 


(Ax-b)T  Q(Ax-b) ,  Q  >  0. 


(6) 


The  design  objective  in  each  case  is  the  minimization  of  the  distance  between 
the  attainable  and  desired  set-point,  where  this  distance  is  defined  by  the 
Euclidean  norm. 


d2  - 


eA  e. 


(7) 


of  the  error  vector  between  x8p  and  the  plant  state  x(t),  i.e., 

e(t)  -  Xgp  -  x(t).  (8) 

An  expression  for  the  dynamics  associated  with  this  error  can  be  derived  by 
differentiating  (8)  and  substituting  in  the  appropriate  terms  from  (2).  The 
result  can  be  expressed  as 


*e(t)  -  Xgp  -  x(t)  -  Ae(t)  -  Bu(t)  -  Axsp  -  Fw(t),  (9) 

where  Axap  represents  the  "set-point  disturbance"  term. 

In  disturbance  accommodating  control  design,  the  control  vector  u(t)  is 
considered  to  be  an  ordered  collection  of  the  various  independent  control 
inputs  which  are  available  to  accomplish  the  primary  control  objective  and  to 
"accommodate"  the  disturbances  which  are  acting  on  the  system.  In  the  design 
of  disturbance  minimization  controllers,  it  is  common  practice  to  split 
(allocate)  the  total  control  u(t)  into  two  parts  as  follows: 

u(t)  -  Up(t)  +  u<j(t),  (10) 

where  up(t)  is  given  the  task  of  accomplishing  the  primary  control  objective 
and  u^(t)  is  given  the  task  of  disturbance  accommodation.  The  part  uj(t)  can 
be  further  subdivided  into  component  vectors,  as  required.  For  the  methods 
considered  in  this  report  u^(t)  will  be  allocated  as 

u^(t)  -  u{j8(t)  +  u<jw(t).  (11) 

The  component  u<j8(t)  will  be  designed  to  accommodate  the  effects  of  the  set- 
point  disturbance  term  while  u,jw(t)  will  be  designed  to  accommodate  the 
effects  of  the  external  disturbance  term.  If  the  plant  is  completely 
controllable  and  is  also  completely  observable,  the  control  up(t)  can  be 
designed  in  the  form 


Up(t)  -  Kx(t) 


(12) 


Given  the  allocation  of  the  control  vector  u(t),  (9)  can  be  re-written 
as 

e(t)  -  Ae(t)  -  BUp(t)  -  Buds(t)  -  Budw(t)  -  Axsp  -  Fw(t).  (13) 

Upon  substitution  from  (4),  (13)  becomes 

kt)  -  Ae(t)-BUp(t)-Buds(t)-Budw(t)-Axsp-FHz(t)-FLxsp-FLe(t).  (14) 

In  the  case  of  (14),  one  would  design  up(t)  in  the  form 

up(t)  — Ke(t)  (15)' 

with  K  chosen  such  that  the  homogeneous  system 

kt)  -  (A+FL+BK) e(t)  (16) 

will  yield  e(t)-*0  "rapidly".  If  one  lets  A  ■  A+FL+BK,  then  (15)  can  be 
expressed  as 

kt)  -  Ae(t)  -  ((A  +  FL)  xsp  +  Buds(t))  -  (FHz(t)  +  Budw(t)).  (17) 

One  of  the  approaches  developed  in  [13]  was  to  provide  minimization  of 
egS  by  use  of  an  allocated  disturbance  control  component.  If  one  assumes  that 
a  unique  steady-state  solution  exists  for  e(t),  i.e.,  all  Xi  of  A  have  nega¬ 
tive  real  parts  and  all  disturbance  terms  have  a  limit  as  time  approaches 
infinity,  then  egS  can  be  found  by  setting *e(t)  in  (17)  to  zero  and  solving 
for  the  e  which  satisfies  the  resulting  equation.  If  this  is  done,  the 
expression  for  the  steady-state  error  is  found  to  be 


egs  *  A-1  t(A  +  xsp  +  Buds  +  FHz«  +  Budw] 
One  now  wishes  to  design  uds  and  udw  such  that 


Sssl 


A-1  [(A  +  FL)xsp  +  Buds] 


and 


£882 


A-1  (FHz  +  Budw) 


(18) 


(19) 


(20) 


are  minimized.  If  Q  in  (6)  is  chosen  to  be  the  identity  matrix  I,  the 
resulting  minimum  norm  control  components  which  minimize  the  Euclidean  norms 
of  (19)  and  (20)  were  shown  ([13])  to  be 


* 

uds 


-<*-i  i/fi 


A  +  FL)xsp 


(21) 
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and 


udw  "  -(A-1-  B  ) ^  A-l  FHz*,  ,  (22) 

where  (.)  denotes  the  Moore-Penrose  generalized  inverse  of  (.).  The  minimal 
norm  set-point  error  egS  was  also  shown  to  be  given  by 


ess  -  U  "  (A-1  B)  (A-1  B)^  [(A  +  FL)xsp  +  FHzJ.  (23) 

Another  approach  developed  in  [13]  was  to  minimize  the  disturbance 
effects  in  (17).  In  this  instance,  ud8(t)  and  u<jw(t)  are  designed  to  minimize 


(A+FL)xsp  +  Buds(t) 


and 


FHz(t)  +  Budw(t) 


(24) 


respectively.  It  was  shown  that,  for  Q=I,  the  minimum  norm  control  components 
which  minimize  the  norms  in  (24)  are  given  by 


uds(t)  -  VU  +  FL)xsp  (25) 

udw(t)  -  -BrFHz( t) .  (26) 

Substitution  of  (25)  and  (26)  into  (17)  results  in  the  following  expression 
for  the  error  dynamics. 


kt)  -  Ae(t)  -  (I  -  BB'SKA  +  FL)xsp  +  FHz(t)].  (27) 

Under  the  same  assumptions  as  to  the  existance  of  a  steady-state  solution  made 
with  reference  to  (17),  the  unique  steady-state  solution  for  (27)  is 

4s  “  A-1  (I  -  BB'*'  )[(A  +  FL)xsp  +  FHzJ.  (28) 

An  example  problem  was  worked-out  in  [13]  using  the  controller  pair 
given  by  (21)  and  (22)  and  the  controller  pair  given  by  (25)  and  (26).  A  com¬ 
parison  of  results  indicated  that  use  of  the  first  controller  pair  did  result 
in  a  smaller  steady-state  error.  In  [14],  additional  data  was  presented  which 
indicated  that,  for  the  same  conditions  considered  in  [13],  the  first 
controller  pair  also  resulted  in  less  transient  excursion  in  e(t). 
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IV 


PLANT  AND  DISTURBANCE  MODELS 


The  plant  state  and  output  models  used  for  the  state  set-point  regulator 
example  of  [13,14]  and  for  examples  to  be  presented  in  this  report  are 


/*i\  u  ll  /xA  fi\  fi\ 

(  I*  {  )+i)u+i#w 

\*2/  L°  W  \2/  W 


y  -  (1,  0)  x 

The  external  disturbance  model  is 


o(t)  , 


i.e.,  the  external  disturbance  was  considered  to  be  constant  between 
a  arrivals.  The  target  state  set-point  vector  is  given  as 


xgp  *  (xspfl»  0*)  •  (33) 

The  plant  given  by  (29)  and  (30)  is  completely  controllable.  For  the 
purpose  of  the  examples,  it  was  assumed  in  [13,14],  and  will  be  assumed  in 
this  report  also,  that  all  necessary  state  information  is  available  from  an 
ideal  reconstructor. 


V 


GEOMETRICAL  CONSIDERATIONS 


Since  the  control  distribution  matrix  B  of  (29)  is  a  2x1  matrix  of  rank 
1,  it  does  not  span  the  state  space,  which  is  two-dimensional  in  this  example. 
Hence,  Bud(t)  will  have  a  limited  set  of  attainable  points  in  the  state  space. 
Also,  the  external  disturbance  distribution  matrix  F  is  of  rank  1  and  thus, 
Fw(t)  will  have  a  limited  range  of  action  in  the  state  space.  As  can  be  seen 
in  Figure  1,  the  lines  of  action  of  Bud  and  Fw  are  not  colinear.  Furthermore, 
as  shown  in  Figure  2,  the  line  of  action  of  the  set-point  disturbance  term,  in 
the  error  state  space,  is  not  colinear  with  the  line  of  action  of  the  control. 
Hence,  no  ud  exists  which  will  completely  absorb  a  non-zero  external  distur¬ 
bance  or  a  set-point  disturbance  resulting  from  a  non-zero  target  state 
set-point. 

Given  that  this  situation  exists  and  that  a  design  objective  is  to  mini¬ 
mize  the  effects  of  the  uncancellable  disturbance,  one  thus  attempts  to  design 
ud  so  as  to  achieve  this  objective  in  some  fashion.  To  illustrate  the  action 
of  the  controllers  shown  in  Section  III,  consider  the  following.  With  respect 
to  the  vectors  Fw^  and  Bu^  shown  on  Figure  1 ,  one  approach  to  the  minimiza¬ 
tion  problem  is  to  first  express  the  vector  Fw  as  the  sum  of  two  component 
vectors,  one  lying  in  the  column  range  space  of  B,  R(B),  and  one  lying  in  the 
orthogonal  complement  to  the  column  range  space  of  B,  R(B)J-.  This  makes  it 
easy  to  see  that  the  component  lying  in  R(B)-^,  which  is  the  component  that  is 
uncancellable,  is  minimized  if  the  component  lying  in  R(B)  is  the  orthogonal 
projection  of  Fw^ ,  onto  R(B).  In  essence,  this  result  is  provided  by  the 
controllers  shown  in  Section  III.  For  instance,  if  udl  is  chosen  as 

udl  —  B^Fwi  (34) 

then 

L  l  (35) 

Budj  +  Fwj  *  -BB  Fwj  +  Fwj  ■  ( I-BB  )Fwj  , 

where  (I-BB^)  is  the  projector  of  Fw^  on  R(B)-*-  along  R(B)  and  BB^  is  the  pro¬ 
jector  of  Fwj  on  R(B)  along  R(B)-*-.  The  uncancellable  part  of  the  external 
disturbance  is  thus  the  component  in  R(B)^-. 


Plant  responses  to  different  disturbance  and  disturbance  control  conditions 
controller  pair  (21)  &  (22) . 
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VII.  EXTERNAL  DISTURBANCE  UTILITY 


A.  Introduction 

In  [13],  the  ’’utility"  U  of  an  external  disturbance  w  was  defined  as 


U 


®ss 


^s 


w  +  0 


(41) 


If  U  is  positive,  then  w  has  aided  in  reducing  the  final  error.  If  U  is 
negative,  then  w  has  further  increased  the  error.  Note  that,  if  the  target 
set-point  is  the  origin,  such  that  the  set-point  disturbance  is  zero,  a  non¬ 
zero,  uncancellable  external  disturbance  will  always  exhibit  a  negative  util¬ 
ity.  It  was  also  shown  in  [13]  that  two  conditions  must  be  satisfied  in  order 
for  an  external  disturbance  to  have  a  positive  utility.  In  order  for  w  to 
provide  a  positive  utility  it  must  satisfy  a  magnitude  condition  given  by 


n  _  _  n 

Z  (a1  +  f1)2  _  E  (ai)2  <  o  (42) 

i-1  i-1 

and  an  angle  condition  given  by 

<2a  +  7,  f>  <  0.,  (43) 


i .  e. , 

90°  <  6  <270°,  (44) 

where 


a  is  the  component  of  Axsp  lying  in  R(B)-*-, 
f  is  the  component  of  Fw  lying  in  R(B)-*-, 

0  is  the  angle  between  the  vectors  2a+f  and  f,  and 
<•>  denotes  the  inner  product. 

It  is  Interesting  to  consider  the  results  obtained  for  the  cases  of 
Figures  3  through  6  when  the  disturbances  are  alternately  set  to  zero  and  the 
disturbance  minimization  control  ud  is  set  to  zero.  Figures  11  to  14  present 
the  results  obtained  when  this  is  done.  Consider  first  the  case  where  ud  is 
designed  to  minimize  £33.  The  separate  effects  due  to  the  external  and  set- 
point  disturbances,  with  no  disturbance  control  active,  for  xap  *  (10.,0.)T 
are  shown  in  Figure  11.  Note  that  when  both  disturbances  are  input,  again 
with  no  disturbance  control  active,  the  resultant  error  is  additive.  In 
Figure  12  however,  where  xsp  *  (-5.,0.)^,  under  the  same  conditions  the  error 
when  both  disturbances  are  present  is  less  than  that  due  to  the  set-point 
disturbance  alone.  In  this  latter  Instance,  the  external  disturbance  has 
aided  in  reducing  the  error,  i.e.,  has  exhibited  a  positive  utility.  Next, 


response  with 


MINIMIZATION  ON 


Figure  8.  Plant  response  with  x  =  (5.,0.)‘  and  ECO)  =  (-10., 0.) 


response  with 


•esponse  with  x  =  (-5.,0.)  and  e(0)  =  (3.,0.) 


response  with  x  =  (-5.,0.)  and  e(0)  =  (5.,0.) 


response 


response 


VI.  PERFORMANCE  COMPARISON  WITH  CONSTANT  DISTURBANCE 

As  was  mentioned  In  Section  III, the  plant  and  disturbance  models  given 
in  Section  IV  were  used  to  work  out  ([13,14])  an  example  problem  in  order  to 
illustrate  the  performance  resulting  from  application  of  the  controller  pairs 
given  by  Equations  (21,  22)  and  (25,  26),  l.e.,  minimizing  the  effects  of  the 
disturbance  on  the  steady-state  error  versus  minimizing  the  effects  of  the 


disturbance  in  the  differential  equation  describing  the  error. 

In  [13],  it  was  shown  that  the  controller  given  by  Equations  (21)  and 
(22)  are 

uds  "  “2*3447xsp,l  (36) 

udw  "  -1.6723m,  (37) 

and  the  controllers  given  by  Equations  (25)  and  (26)  are 

Ud8  ■  ”®*2xsp,l  (38) 

* 

udw  “  “0.6*.  (39) 

The  stabilization  control  up,  designed  according  to  Equations  (15)  and  (16), 
was  chosen  to  be 

up(t)  -  -Ke(t)  -  (3.,  0.36)e(t).  (40) 


In  [14],  results  were  presented  which  show  the  transient  behavior  of  the 
error  for  each  pair  of  controllers  for  a  case  with  e(0)  ■  (10.,0.)T.  Figure  3 
presents  these  results.  As  can  be  seen  from  Figure  3,  for  the  conditions 
assumed,  designing  uj  to  mlmimize  the  steady-state  error  resulted  in  the 
smaller  steady-state  error  and  also  in  a  better  transient  excursion 
performance. 

It  is  of  interest  to  know  if  the  controllers  designed  as  in  Equations 
(36)  and  (37)  will  always  result  in  the  best  overall  performance.  This 
controller  design  technique  will  always  result  in  the  smallest  steady-state 
error  which  it  is  possible  to  achieve  for  a  given  set  of  conditions,  l.e,  a 
given  target  set-point.  Initial  conditions  and  constant  external  disturbance, 
since  it  results  in  controllers  which  are  designed  specifically  to  minimize 
the  steady-state  error.  But  what  about  the  transient  excursion  performance? 

In  order  to  examine  this  question,  several  additional  cases  have  been 
simulated.  Results  are  presented  in  Figures  4  to  10.  In  all  cases,  the 
controller  pair  designed  to  minimize  egS  did  give  the  smallest  steady-state 
error.  Neither  controller  pair  exhibited  a  consistent  advantage  in  transient 
excursion  performance;  however,  the  performance  associated  with  controller 
pair  (36,  37)  was,  in  all  cases,  at  least  as  good  as  that  associated  with 
controller  pair  (38,  39). 


action  for  disturbances  and  control 


action  of  external  disturbance  and  control 


consider  Che  case  where  u<j  is  designed  to  minimize  Che  disturbance  effects  on 
e.  The  results  for  this  case  are  shown  in  Figures  13  and  14.  Again,  it  can 
be  seen  that  for  x8p  ■  (-5.,0.)T  the  external  disturbance  assists  in  reducing 
the  total  error. 


This  Section  contains  two  examples  which  show  the  region  of  positive 
utility  of  external  disturbances,  the  geometry  Involved  in  the  state  space  and 
several  other  Interesting  results  for  the  set-point  regulator  model  of  Section 
IV*  The  first  example  is  for  the  case  with  x8p  -  (10.,0.)T  and  the  second  is 
for  the  case  with  x8p  -  (-5.,0.)T. 

B.  Example  1 


For  this  first  example,  the  target 
(10.,0.)T.  The  set-point  disturbance  vector 


state  set-point  is  x, 
is  thus  given  by 


(45) 


The  external  disturbance  vector  in  state  space  is  given  by 


(46) 


Figure  15  is  a  plot  of  the  state  space  showing  the  two  disturbance  vectors  and 
the  line  of  action,  R(B) ,  of  the  control. 


As  was  mentioned  in  Section  V,  each  of  the  disturbance  vectors  can  be 
expressed  as  the  sum  of  two  component  vectors  [13],  one  lying  in  R(B)  and  one 
lying  in  R(B)^-.  The  component  of  each  vector  lying  in  R(B)  can  be  found  by 
multiplying  the  vector  by  BB^  and  the  component  lying  in  R(B)-*-  can  be  found  by 
multiplying  the  vector  by  (I-BB^).  The  components  in  R(B)-^  are  then  used  in 
checking  for  satisfaction  of  conditions  (42)  and  (44).  Instead  of  (42),  one 
can  also  use  the  Inequality  [13], 


f 


< 


(47) 


to  establish  an  upper  bound  on  the  magnitude  which  f  can  have  and  still  allow 
w  to  exhibit  a  positive  utility. 


For  the  set-point  disturbance  of  (45),  the  component  lying  in  R(B)-*-  is 
found  to  be 


a  -  (I  -  BB  ;  Ax 


sp 


0.8  -0.4 

-0.4  0.2 


and  the  component  of  (46)  lying  in  R(B)^  is 

j  f  0.8  -0.4 

f  -  (I  -  BBr)FHz  -  I 


/l0.\  /  8. 

V  0.1  \-4. 

(I)  •  C:) 


(48) 


-0.4 


0.2 


(49) 


ipnses  to  different  disturbance  and  disturbance  control  conditions 
pair  (25)  &  (26) . 


These  components  are  also  shown  In  Figure  IS.  As  can  be  seen,  they  act  in  the 
same  direction  and  hence  will  reinforce  each  other,  l.e.,  one  would  expect  w  to 
exhibit  a  negative  utility  in  this  case.  Substituting  from  (48)  and  (49)  into 
(42),  one  obtains 


2  _ 


125.  -  80.  >  0. 


(50) 


and  it  is  evident  that  9  ■  0*;  therefore,  w  does  not  satisfy  the  conditions 
for  0>0. 


According  to  the  requirement  given  by  (44),  for  a  positive  utility 
in  this  example  7  must  lie  to  the  "left"  of  R(B).  Since  the  line  of  action  of 
w  is  a  line  with  a  slope  of  +1  in  the  state  space,  those  external  disturbances 
which  may  exhibit  a  positive  utility  must  be  negative.  Furthermore,  the  bound 
on  the  allowable  magnitude  of  7,  as  indicated  by  (47),  specified  that 


2a 


17.89 


Since  f  -  (I-BB  )FHz,  one  has  that 


(51) 


(52) 


(53) 


z 


< 


40. 


(54) 


For  w  to  exhibit  a  positive  utility  when  xSp  -  (10., 0.)T,  FHz  must 
be  negative,  with  the  magnitude  of  z  less  than  40.  Figure  16  indicates  the 
regions  of  positive  and  negative utility  for  possible  external  disturbance 
magnitudes.  Since  FHz  is  a  line  in  state-space,  the  utility  regions  of  w  are 
as  follows: 


w  >  0. 
w  -  0. 
-40.  <  w  <  0. 

w  <  -40. 


U  <  0. 
U  -  0. 
U  >  0. 
U  <  0. 


(55) 


Figure  17  shows  eS8  as  a  function  of  external  disturbance,  with  and 
without  disturbance  minimization  control  and  with  ud  -  0.  Figure  18  is  an 
expanded  scale  plot  of  the  region  in  Figure  17  from  w  -  -20.  to  w  -  -10. 
Figure  19  shows  the  percent  reduction  in  e8S  achieved  by  the  two  disturbance 
minimization  controllers  and  Figure  20  is  an  expanded  plot  of  a  portion  of 
Figure  19  for  w  between  -20.  and  -10.  Note  from  Figure  18  that  there  is  a 
range  of  values  for  w  for  which  the  controller  designed  to  minimize  distur¬ 
bance  effects  on  e  causes  an  increase  in  e8s  from  what  would  be  obtained  if 
ud  =0.  Also,  note  that  there  is  one  value  for  w,  with  ud“0.,  which  equals 
the  performance  of  the  controller  designed  to  minimize  egS* 
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Geometry  of  disturbances  relative  to  the  line  of  action  of  the  control 


external  disturbance  magnitude 


VW*  WTWm  H  L  ■V-’.'Tk  C'VTTFTTE  T<.T»."V  T  V  \ -.’  W  ’  '  mV'  \  «.v  VT^i  '  V  V  '  ■_'  '/  *JT 


C.  Example  2 


In  this  example,  the  target  state  set-point  is  xgp  «  (-5.,0.)T.  The 
set-point  vector  is  thus  given  by 


Ax 


sp 


1. 


(56) 


The  external  disturbance  vector  is  as  shown  in  (46).  The  component  of 
Axsp  lying  in  R(B)-*-  is  found  to  be 


0.8  -0.4 

-0.4  0.2 


(57) 


and  the  component  of  Fw  in  R(B)-l  is  as  given  by  (49).  The  geometry  of  the 
disturbance  vectors  and  the  components  in  R(B)1  for  this  case  are  shown  in 
Figure  21. 


As  can  be  seen  from  Figure  21,  the  uncancellable  components  of  the  set- 
point  and  external  disturbances  act  in  opposite  directions,  unlike  the 
situation  which  existed  in  the  first  example,  and  one  would  expect  w  to  exhi¬ 
bit  a  positive  utility.  Upon  substituting  from  (49)  and  (47)  into  (52),  one 
obtains 


a  +  f 


2 


5. 


20. 


-15.  <  0. 


(58) 


and  it  is  evident  from  Figure  21  that  0  -  180*;  therefore,  w  does  satisfy  the 
conditions  for  positive  utility. 


The  bound  on  the  all'’  >le  magnitude  of  f  for  positive  utility  is 


2a 


8.944, 


(59) 


so. 


20. 


(60) 


For  U  >  0.,  with  xgp  -  (-5.,0.)T,  one  must  have  9Oo<0<27O,,  and  ||  z  ||  <20. 
Figure  22  shows  the  region  of  positive  utility  in  the  state  space.  Figure  23 
shows  the  error,  as  a  function  of  external  disturbance,  with  and  without 
disturbance  minimization  control.  Figure  24  shows  an  expanded  scale  plot  of 
the  portion  of  Figure  23  between  w  ■  0.  and  w  *  10.  Figure  25  shows  the  per¬ 
cent  reduction  in  Cg8  obtained  by  use  of  the  disturbance  minimization 
controllers  over  the  case  when  u<j-0. ,  and  Figure  26  is  an  expanded  scale  plot 
of  the  portion  of  Figure  25  between  w  ■  0.  and  w  -  10. 


As  can  be  seen  from  Figure  24  there  is  again  a  region  for  w  in 
which  the  controllers  designed  to  minimize  disturbance  effects  on  e  cause 
larger  eg8  than  would  Uj=0.  There  is  also  a  value  of  w,  with  u^a0.,  which 
equals  the  performance  of  the  controllers  designed  to  give  the  minimum  eaa. 
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0001 

COMMON  X(14). DX( 14).  KUTT4.  DT,  NX 

0002 

DIMENSION  XDAT ( 50 > 

0003 

NVAR-6 

0004 

UR  I TE ( 20 ) NV AR 

0005 

W— 42. 

0006 

DO  5050  IJ-1. 120 

0007 

W-W«-l. 

OOOS 

DO  100  1-1.  14 

0009 

DX(I>-0. 

0010 

X(I>»0. 

0011 

100  CONTINUE 

0012 

TIME-0. 

0013 

XSP— 5. 

0014 

NX-16 

0015 

Xl-5. 

00l4> 

X3-5. 

0017 

c 

*****  INITIAL  CONDITIONS  ON  PLANT  STATES  ***** 

ooia 

X(1)-X1 

0019 

X (3J-X1 

0020 

X ( 5>— XI 

0021 

X(7)» XI 

0022 

X(9)»X3 

0023 

X ( 11 )— X3 

0024 

X ( 13>-X3 

0025 

C 

******************************-»•»************* 

0026 

DT— 0.  05 

0027 

UDS1S-0. 

00 28 

UDW1S-0. 

0029 

UDS2S-0 

0030 

UDW2S-0. 

0031 

UD1T-0. 

0032 

UD2T-0. 

0033 

IPRT-0 

0034 

1000  CONTINUE 

0035 

IF  (TIME.  GE.  10.  >  00  TO  9999 

0036 

IPRT-IPRT+1 

0037 

DO  200  KUTTA-1.4 

0038 

C 

*****  DIST  MINIMIZING  CONTROL  VECTORS  FOR  MIN  NORM  OF 

***** 

0039 

/*» 

w 

*****  STEADY-STATE  ERROR 

***** 

0040 

UDS1— 2.  3447*XSP 

0041 

UDU1  —  1.  6723*W 

0042 

c 

****************************-»  ************************* 

0043 

C 

*****  DIST  MINIMIZING  CONTROL  VECTORS  FOR  MIN  DISTUR- 

***** 

0044 

r 

*****  3ANCE  CONTRIBUTION  IN  CI--E?ENTIAL  EG  FOR  ERROR 

***** 

0045 

UDS2— 0.  2*XSP 

0046 

UDU2--0.  6*W 

0047 

c 

************************************************************** 

0048 

c 

*****  PLANT  DIFF  EOS  WITH  UD  T"  MIN  NORM  OF  STEADY- 

***** 

0049 

c 

*****  STATE  ERROR 

***** 

0050 

DX(  1  >— 2.  *X(  1  >+0.  64*X  (  2>-XSP— UDS1-W— UDW1 

0051 

DX (2>— 6.  *X(1  >*0.  28*X(2)-2.  *UDSl-W-2.  *UDW1 

0052 

c 

*******************************i-»************************** 

0053 

c 

*****  PLANT  DIFF  EQS  WITH  UD-0 

***** 

0054 

DX  (3)— 2.  *X  ( 3 )  *0.  64*X  ( 4  > -  ‘ SP— W 

0055 

DX ( 4)— 6.  *X (3)  +0.  28*X(4)-w 

0056 

c 

*****************************  »•>•>■*  ***************•:♦********** 

0057 

c 

*****  PLANT  DIFF  EOS  WITH  UDW-* 

***** 

A-2 


APPENDIX 

DIGITAL  SIMULATION 
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IX 


SUMMARY 


It  was  shown  in  Section  VI  that  u<j  designed  to  minimize  steady-state 
error  does  result  in  a  smaller  CgS  than  does  u<j  designed  to  minimize  distur¬ 
bance  effects  on  e.  However,  neither  design  exhibited  a  consistent  advantage 
in  reducing  transient  excursions  in  the  plant  response. 

Two  examples  were  worked  out  in  Section  VII  in  order  to  illustrate  the 
geometry  associated  with  disturbance  minimization  for  a  second  order  state 
set-point  regulation  problem  and  the  regions  of  positive  and  negative  utility. 
These  examples  demonstrate  the  fact  that  it  is  possible  for  an  uncancellable 
external  disturbance  to  assist  in  reducing  the  set-point  error.  From  Figures 
17  and  23,  it  can  be  seen  that  this  is  true  whether  or  not  the  disturbance 
minimizing  control  is  present.  The  disturbance  minimizing  control  did  provide 
a  smaller  achievable  steady-state  error  and  did  extend  the  range  of  external 
disturbances  for  which  a  positive  utility  could  be  obtained.  This  can  be  seen 
as  illustrated  in  Figure  17.  For  the  case  with  u^  ■  0.  the  steady-state  error 
for  -30.<w<0.  was  less  than  or  equal  to  the  steady-state  error  for  w  *  0. 
However,  for  the  case  where  u<j  was  designed  to  minimize  e8S,  the  range  of 
values  of  w  for  which  U>0  goes  to  -40.  It  is  also  apparent  from  Figures  17 
and  23  that  application  of  the  disturbance  minimization  controller  can  result 
in  zero  steady-state  error,  for  a  particular  value  of  w  in  each  case,  i.e.,  w 
■-2Xgp  j.  This  result  can  be  verified  from  manipulation  of  Equation  (23). 

In  Section  VIII,  it  was  shown  that  it  is  possible  to  reduce  control 
energy  expenditure  in  cases  where  w  provides  a  steady-state  error  equal  to 
that  given  by  u<j  designed  to  reduce  egg. 


Figure  31.  Control  energy  vs.  time,  corresponds  to  case  in  Figure  30. 


■MM 


error  for  various 


To  illustrate  use  of  this  property.  Example  1  was  simulated  with  the 

magnitude  of  w  changing  at  three  second  intervals  up  to  a  time  of  9  seconds. 

Figure  27  shows  the  set-point  error  obtained  when  no  disturbance  minimization 
control  is  applied.  Figure  28  shows  the  set-point  error  obtained  when  u^  * 
u^  except  during  the  interval  when  w*-14.,  at  which  time  u^  is  set  to  zero. 
The  control  energy  corresponding  to  this  case  is  shown  in  Figure  29.  Figure 
30  shows  the  set-point  error  when  u<j  *  uj  for  all  values  of  w  and  Figure  31 
shows  the  corresponding  control  energy.  A  comparison  of  Figures  28  and  30 
shows  that  the  two  curves  are  identical.  A  comparison  of  Figures  29  and  31 

shows  that  the  control  energy  expenditure  is  reduced  by  about  140  units  when 

U4  is  set  to  zero  for  w— 14.  Of  course,  these  results  were  obtained  under 
the  assumption  that  w  is  known  exactly  and  that  the  time  of  occurence  of  each 
change  in  w  is  known  exactly,  but  they  serve  to  illustrate  the  possibilities. 


VIII.  REDUCING  CONTROL  ENERGY  EXPENDITURE 


As  can  be  seen  in  Figures  18  and  24,  there  is  one  point  where  the 
steady-state  set-point  error  is  the  same  without  u^  as  it  is  with  u<j  designed 
to  minimize  ega.  That  is,  for  a  given  target  set-point  there  is  an  external 
disturbance  w  for  which  the  controller  u^  could  be  set  to  zero  without  any 
performance  degradation.  There  is  also  a  small  interval  about  this  particular 
w  within  which  w  could  be  located  and  for  which  u<j  would  contribute  little  to 
further  error  reduction.  This  would  permit  u<j  to  be  set  to  zero,  in  order  to 
conserve  control  energy,  without  deteriraental  effects  on  performance.  Note 
that  the  particular  w  is  close  to,  but  not  coincident  with,  the  point  of  mini¬ 
mum  Ggg  for  the  u<j*0.  curve  in  each  case. 

In  order  to  determine  the  relationship  between  w  and  xgp  at  these  equal 
performance  points  for  the  plant  and  disturbance  models  of  Section  VII,  one 
would  equate  the  norms  of  Equation  (23)  and  the  equation  resulting  from  (18) 


when  U(js  *  u,jw  *  0.,  i.e.,  the  norms  of 

4a  «  [I  -  (A“1B)(A“1B)’6  ]A-1  (Axsp  +  Fw)  (23) 

and 

EgS  -  A-l  (Axsp  +  Fw).  (61) 

The  norm  squared  of  (23)  is  given  by 

||4sj|2  *  °*8xsp,l  +  0.8xsptl  w  +  0.2w2  (62) 

and  the  norm  squared  of  (61)  by 

Egg  2  -  3.352xgp>1  +  4.441xsp>1w  +  1.498w2.  (63) 

If  one  equates  (62)  and  (63),  the  resulting  expression  is 

x2p,i  +  1.426Xgp>1w  +  0.508w2  ■  0  (64) 

and  the  solution  for  w  is  given  by 

w  -  -1.402xsp>|.  (65) 


For  a  given  target  set-point  of  the  form  xsp  -  (xsp  0.)T,  Equation  (65) 
will  give  that  value  for  w  which,  if  it  occurs,  will  enable  the  control  uj  to 
be  set  to  zero  without  loss  of  performance.  In  Example  1,  the  value  of  w 
from  (65)  is  w  -  -14.  In  Example  2,  the  value  of  w  from  (65)  is  w  -  -7.0.  It 
can  be  seen  from  Figures  18  and  24  that  these  are  the  correct  values. 


XTERNRL 


08 '2 


i - r 

8  08  '9 


08‘t 


(3ms-AQV3iS)  3amiN9vw  aoaa3 


0058  DX<5>  —  2.  *X<5/+0.  64*X<6>-XSP-W-UDS1 

0059  DX<6>— 6.  *X(5>+0.  28*X(6)-w-2.  *UDS1 

0080  C  *****************♦********»+*•»•■»-►-*•■»■*■■»■***■»***■**■***■*■**•***♦**** 

0081  C  *****  PLANT  DIFF  EOS  WITH  UDS-C.  ***** 

0082  DX<7>  — 2.  *X<7>+0.  64*X(8>-XSP-W-UDW1 

0083  DX<8>— 6.  *X(7)+0— 0.  28*X(S>— W— 2.  *UDU1 

0084  C  »*»***»»*»**»*»»»***»»»»»»»**»****(»*»*»*»»»»»»**»»*»»*»»»»» 

0085  C  *****  PLANT  DIFF  EOS  WITH  UD  TO  MIN  DIST  CONTRIBUTION  ***** 

0068  C  *****  IN  DIFF  EOS  FOR  ERROR  ***** 

0067  DX<9>—  2.  *X (9)  +0.  65*X ( 10  > -XSP-UDS2-W-UDW2 

0068  DX  <  10  >  —6.  *X  ( 9 )  *0.  28«X  ( 10 )  -W-2.  *UDS2-2.  *UDW2 

0070  C  *****  PLANT  DIFF  EOS  WITH  UDW-O.  ***** 

0071  DX(ll)— 2.  *X<  13)  +0.  64*X<lii-XSP-W-UDS2 

0072  DX(12>— 6.  *X(13)*0.  28* X < l 4 > -W-2.  *UDS2 

0073  C  ******************************  .>■•»•>-»»*»»*»♦»»»»♦»»»*»»»»»»»»»♦ 

0074  C  *****  PLANT  DIFF  EOS  WITH  UDS-C  ***** 

0075  DX< 13)*-2. *X(15>+0.  64*X i la >-XSP-W-UDW2 

0076  DX ( 14)— 6.  *X(l5)+0.  28* X ( 16) -W-2.  *UDW2 

0077  C  ****************** ■**♦♦*******■» -*-»-t************************** 

0078  00  TO  (30. 80. 30.  40) . KUTTA 

0079  30  CONTINUE 

0080  TIME-TIME*.  5*DT 

0081  40  CONTINUE 

0082  60  CALL  RUNK 

0083  200  CONTINUE 

0084  C  *****  MEASURE  OF  CONTROL  ENERGV  ***** 

0085  UDS1S— UDS1S*ABS ( UDS1 »DT ) 

0086  UDW1S— UDW1S*A3S(UDW1*DT) 

0087  UDS2S— UDS2S*AB  S ( UDS2*DT ) 

0088  UGW2S— UDW2S*ABS ( UDW2*DT ) 

0089  UD 1 T -UDS 1 S*UDW 1 S 

0090  UD2T -UDS2S+UDW2S 

0091  C  ************************************************************* 

0092  C  *****  ERROR  MAGNITUDE  IN  EACH  CASE  ***** 

0093  ERR 1-SQRT ( X ( 1) **2*X ( 2 ) **2 ( 

0094  ERR2-SQRT ( X ( 3 ) **2*X ( 4 ) **2 ) 

0095  ERR3-S«T(X(5)**2*X(6)**2; 

0096  ERR4-SQRT ( X ( 7 ) **2*X ( 8  >  ** 2 : 

0097  ERR5-9QRT  ( X  ( 9 )  **2*X  ( 1 0  >  *<-2 ) 

00«»S  ERR6-SGRT  ( X  <  1 1  >  **2*X  (12;  -'■•'2 > 

0099  ERR7-SGRT ( X ( 13 ) »*2*X (14/ --2 > 

0100  0  ***************************************************************** 


0101 

IF < IPRT.  NE  2i  00  TO  500 

0102 

IFRT-0 

0103 

500 

00  TO  1000 

0104 

50 

FORMAT  ('  TIME  F8.  4.-  ERP.l 

FS.  4,  ' 

0103 

. 

'  ERR3  F3.  4.  '  ERR 4  ’ .  F5.  4. 

'  ERR5  ' 

01  Co 

9999 

CONTINUE 

0107 

c  ***** 

PERCENT  REDUCTION  FROM  CASE  WITH 

UD-0 

♦**r** 

0108 

XPCES- <  <  ERR2-ERR 1 ) /ERR2 > *1 00. 

0109 

XPCED«( (ERR2-ERR5) /ERR  2; -  CC. 

0110 

C  *********  #*->**************-»—*•. . * . ►***-v**-!’**************** 

'111 

c  ***** 

PLOT  VARIABLES 

. 

***** 

Ol  i2 

XDAT ( 1 ) »W 

0113 

XDAT (2  >«XPCE3 

Oil* 

XDAT (3)-XPCED 

on  3 

XDAT ( 4) *ERR 1 

0U6 

XDAT i 3) -ERRS 

0il7 

XDAT(6)«ERR3 

01  is 

C  #»»♦»»»««»»«»»»»#«♦♦♦«»»♦»>*■»•“■•■»-*■»«»•»♦♦»♦«»♦«»**♦♦**«**«» 

0119 

UR  ITE(  20 )  (XDATCI).  I-I.N.aR) 

GliO 

PRINT  *.  IJ.  W 

0121 

3030 

CONTINUE  • 

01 22 

9T0P 

a 

0123 

END 

PROGRAM 

SECTIONS 

- 

Name 


Attribute* 


0  4C0DE 

ills 

PIC 

CON  REL 

LCL 

SHR  EXE 

RO 

2  *LOCAL 

324 

PIC 

CON  REL 

LCL 

NOSHR  NOEXC 

RO 

3  * BLANK 

124 

PIC 

OVR  REL 

OBL 

SHR  NOEXE 

RO 

Total  Space  Allocated 

1566 

i TR r  POINTS 

Address 

Type 

Name 

0-0000:000 

DMPLIS4MAIN 

-R IABLES 

Address 

Type 

Naina 

Aadress 

Type 

Nam* 

Address 

Tyi 

3-00000074 

R*4 

DT 

2-0000-114 

R*4 

ERR1 

2-00000118 

R 

2-00000120 

R*4 

ERR4 

2-00000124 

R*4 

ERRS 

2-00000 12S 

R- 

2-00000004 

1*4 

I 

2-00C0C0D0 

1*4 

IJ 

2-00000100 

I< 

2-Q00000CB 

1*4 

NVAR 

3-00000078 

1*4 

NX 

2-000000DB 

R- 

2-0000CCFC 

R*4 

UD2T 

2-00000104 

R*4 

UDS1 

2-COOOOOES 

R’ 

2-OCOOOOFO 

R*4 

UDS2S 

2-0000C 109 

R*4 

UDW1 

2-000000EC 

R* 

2-0C0000F4 

R*4 

U0U2S 

2-0C00CCCC 

R*4 

U 

2-000000E0 

R- 

2-00000134 

Re* 

XPCED 

2-coo: 1120 

R*4 

XPCES 

2-OOOOOODC 

P- 

ARRAYS 


Address 

Type 

Nam* 

3^  res 

Dimensions 

3-00000038 

R*4 

OX 

36 

(14) 

3-00000000 

R*4 

X 

36 

<  14  > 

• 

2-OOOOCOOO 

R*4 

XDAT 

ZOO 

(30) 

A-4 


/crx^yjC'^y^t,^ 


UK 


SUBROUTINE  RUNK 

COMMON  X ( 14) .  DX  < 14 ) ,  KUTTA.  DT ,  NX 
DIMENSION  XA( 14),  DXA< 14> 

00  TO  <10. 30- SO. 70).  KUTTA 
DO  20  1*1.  NX 
XAC I )*X ( I ) 

DXA ( I ) *DT *DX  < I > 

X< I »— X< I >♦.  3*0XA( I ) 

RETURN 

TDT-2.  *DT 

HOT*.  S*DT 

DO  40  1*1,  NX 

DXA< I )-DXA< I )+TDT*DX  < I ) 

X ( I >«XA< I )+HDT+DX < I > 

RETURN 
DO  60  1*1.  NX 
VDT*DT*DX < I ) 

DXA< I ) *DXA( I ) +2.  *VDT 
X  < I ) *XA  < I ) *VDT 
RETURN 
DO  80  I-l.NX 

X( I )*XA< I )*(DXA< I )+DT*DX ' I >  > /6. 

RETURN 

END 


rSjORAM  SECTIONS 


Attr ibut*s 


0  *CODE 

2  SLOCAL 

3  *3 LANK 


PIC  CON  REL  LCL  SHR  EXE 
PIC  CON  REL  LCL  NOSHR  NOEXE 
PIC  OUR  REL  OBL  SHR  NOEXE 


Total  Space  Allocatad 


ENTRY  POINTS 

Address  Typ*  Nam* 


3- JOOCOCOO 


variables 


Address 

Type 

Nam* 

Add?  ss  s 

Typ* 

Nam* 

Address 

Ty 

3-00000074 

R*4 

DT 

2-000 0007S 

R*4 

HDT 

2-00000070 

I 

3-0000QC78 

1*4 

NX 

2-00000074 

R*4 

TDT 

2-0000007C 

F 

DISTRIBUTION 


Copies 


US  Army  Materiel  System  Analysis  Activity 
ATTN:  AMXSY-MP 

Aberdeen  Proving  Ground,  MD  21005  1 

Commander 

US  Army  Research  Office 
ATTN:  AMXROPH,  Dr.  R.  Lontz 
P.  0.  Box  12211 


Research  Triangle  Park,  NC  27709  1 

Dr.  C.  D.  Johnson 

4001  Granada  Drive 

Huntsville,  AL  35802  .  1 

AMSMI-R,  Dr.  McCorkle  1 

-R,  Dr.  Rhoades  1 

-RG,  Dr.  Yates  1 

-RGN,  Mr.  McLean  10 

-RGS,  Dr.  Kelly  1 

-RGT,  Mr.  Jordan  1 

-RPR  15 

-RPT  (Record  Copy)  1 

-LP,  Mr.  Voigt  1 


END 

FILMED 

9-85 

DTIC 

!SN3cX?  V’V’!Ar  L  !%*  rl3  'prO*H3kf 


tr"  •  v*  •«>'.  -?  >  '*:•  r. 


